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We have computed in [hep-th/1 107.4320] the glueball spectram in a certain sector of the large-A' 
YM theory by solving by a change of variables the holomorphic loop equation for cusped twistor 
Wilson loops supported on certain Lagrangian submanifolds and by evaluating the correlators 
of surface operators supported on these Lagrangian submanifolds. We have shown that the cor- 
relators of composite surface operators of length L reproduce in the large-L limit the leading 
logarithms of perturbation theory of the corresponding glueball propagators, including the correct 
anomalous dimensions. In this paper we show that the correlators of surface operators match in 
the large-L limit the stronger constraints arising by the operator product expansion, according 
to Migdal technique of computing the spectral sum over the glueball including the subleading 
asymptotics given by the Euler formula. Finally, we discuss correlators of surface operators for 
finite L. 
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1. Introduction 

The aim of this paper is to show that the recently computed glueball propagators [|l|] in a 
certain sector of the \wge.-N YM theory reproduce not only the logarithms of perturbation theory, 
including the anomalous dimensions, but satisfy also the more stringent Migdal criterium ^ about 
the operator product expansion {OPE). 

In fact, more precisely, our statement refers to an asymptotic identification for large L, dis- 
cussed in ^ and recalled in the following section, between correlators of composite surface opera- 
tors of length L supported on certain Lagrangian submanifolds and the corresponding correlators of 
local composite single trace gauge invariant operators with momentum dual, in the Fourier sense, 
to the Lagrangian submanifold. 

At the end of next section we discuss also correlators of surface operators for finite L. 

Another point of this paper is to discuss the results ^ for the glueball propagators in the light 
of the philosophy of Migdal "meromorphization" 

In this paper we limit ourselves to analyzing the final results without recalling their computa- 
tion, which is extensively explained in 

But we recall in this introduction some basic conjectures and results of [^, ^. 

The first basic conjecture by Migdal, that dates long ago [^, is that connected two-point 
correlation functions of local single trace gauge invariant operators, 0{x), of \wge.-N YM or QCD 
are saturated by a sum of pure poles (for simplicity of notation we consider the scalar case only): 



/ 



e'P' < O{x)O{0) >,onn d'x = Go{p^) = £ (1-1) 



Zk{p^) are supposed to be analytic functions of with no poles, positive at the poles of Eq.(l.l), 
after analytic continuation to Minkowski. 

The conjecture is based on the estimate that the interaction is suppressed by powers of ^ in the 
expansion by planar graphs that defines the \wge.-N limit of YM or QCD [Q] and on the assumption 
that the theory confines the chromoelectric charge, is such a way that the physical spectrum at the 
leading large-A'^ order is made by infinite towers of free particles, the glueball and the mesons. 

In the string picture these particles would arise as the modes of a fluctuating string, open for 
the mesons and closed for the glueball (see ^ and references therein). 

Naively the two point function, Go{p^), has to satisfy the renormalization group {RG) equa- 
tion: 

{-^+m^+^7o{8))Go{p^) = (1.2) 
dlogp dg 



where yo {§) is the anomalous dimension 



and j8 (g) the beta function 



7»te) = ^ (13) 



dlogp 
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More recently it has been argued in [Q] (v.l) that only the spectral function (Minkowski signature 
is understood): 

Y^Z,{ml)8{p^-ml) (1.5) 

k 

has in fact to satisfy this property, because of the existence of additive renormalizations of the two- 
point function of composite operators and according to the general idea that only the locations of 
the poles and their residues have a physical meaning. Therefore it is not restrictive to write [Q] : 

Go{p') = l^^^, (1.6) 

This is the second basic result of Migdal "meromorphization" [^. 

It is quite clear [^] why the tower of mesons and glueball has to be infinite, in order to match 
the perturbative computation of the correlators. 

Indeed for an operator, 0{x), of naive dimension L ^: 

Go{p^)-Zl{p^)p^^-'\og{^) (1.7) 

in perturbation theory. Thus the sum in Eq.(1.6) can reproduce, asymptotically for large p^, 
Eq.(1.7) only if it is infinite. 

In perturbation theory the power of momentum, p^^^'^, arises by naive dimensional analysis. 
The log(^) factor would imply conformal behavior of the correlator at lowest order. The factor of 
Zq{p^) accounts for the anomalous dimension of the operator 0{x). Zo{p^) is given at one-loop 
order in perturbation theory by 

2 

Zo(/) = l + |/log(^) (1.8) 

for some pure number y. The anomalous dimensions are scheme independent only at one loop. 
The one-loop RG improved expression for Zo{p^) is: 

Zo(/)~[log(^)]* (1.9) 

Recently ^ it has been pointed out in |Q] that the spectral sum occurring in the glueball propagator 
Go{p^) has to satisfy, in addition to Eq(1.7), stronger constraints that arise by the OPE. 

The Migdal criterium for the glueball propagators can be described as follows. 

If the discrete spectral sum on the glueballs is known, it can be evaluated by the Euler formula 



Go{p') =Y^Gu{p') 
k 

B 



poo ^ D 

G,{p')dk-l^^ji{d,y-'[G,{p')]^^^ (1.10) 



'The symbol ~ stays for "equal up to constant irrelevant numerical factors" or "equal up to irrelevant additive terms" 
depending on the framework. 

^We understand that Migdal argument has been known to him and to his collaborators for long, but never published 
before Jl. 
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where Bj are the Bernoulh numbers. 

The leading contribution to the spectral sum arises substituting to the discrete sum the integral. 
This has to reproduce the perturbative result, Eq.(1.7). 

However, Euler formula furnishes the subleading asymtotics to the integral, that contain in- 
verse powers of the momentum p^, because of the derivatives. These inverse powers are interpreted 
as non-perturbative contributions arising by condensates of higher dimension operators in the OPE 
of the Fourier transform of O{x)O{0). 

In iQ] it is remarked that also these terms in general have to carry anomalous dimensions, 
otherwise these terms would imply the existence of unknown conserved currents in YM or in QCD, 
in addition to the traceless part of the stress energy tensor. 

In this paper we check our results [|]] for the glueball propagators against Migdal criterium. 

However, we should mention that there is an alternative school of thought ^ that does not give 
any physical meaning to the condensates of higher dimension operators, because of the ambiguity 
associated to additive renormalizations, necessary to make sense of them in perturbation theory. 
According to this school of thought the only constraint that can be meaningfully tested in pertur- 
bation theory is the leading perturbative term given by Eq.(1.7), the condensates being to every 
order of perturbation theory, because they are proportional to powers of the RG invariant scale: 

Aecz) = Aexp(-— i^)(j8o4cz))~^(l + -) d-H) 

^mQCD 

the OPE being ambiguous by the aforementioned additive renormalizations and by non-perturbative 
terms vanishing in perturbation theory. 



2. glueball propagators 

Recently we computed in [|l|] the glueball spectrum in a certain sector of the large N YM theory, 
by reducing to a critical equation a new kind of loop equation of \wge.-N YM for special cusped 
Wilson loops, called cusped twistor Wilson loops. 

Twistor Wilson loops are supported on Lagrangian submanifolds of space-time and the corre- 
lators that can be computed solving the new loop equation are built by extended objects supported 
on these Lagrangian submanifolds, known as surface operators. 

The crucial property that allows solving the new loop equation for a twistor Wilson loop, as 
opposed to an ordinary Wilson loop, is the absence of cusp anomaly in the \wge.-N limit for the 
former. This property plays a role in the interpretation of the correlators of surface operators. 

In fact a whole family of correlators of the Fourier transform of composite operators of naive 
dimension 4L, 0^{p^,p-), constructed by surface operators supported on the aforementioned La- 
grangian submanifolds, are computed in [|l|] : 

£(A)8i-8A,'2^2p4 f <trMTr^,{^fl)'^{x+,x^)trf,Tr^,{^flf (0,0) J^x 
~< TryrO^{p+,p-)Tr_yO^{-p+,-p-) >fln 

(2.1) 

^We would like to thank Massimo Testa for pointing out to us this different point of view. 
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For large L there is an identification [|l|] between the Fourier transform of our composite surface 
operators and the Fourier transform of composite local operators 

0^{p+,p-) ~ \{F^,+iF^;)\''^{p+,p^,p+,p^) {1.1) 

based on the coincidence of quantum numbers and anomalous dimensions at large L described 
below. 

^ap ^^'^ anti-selfdual {ASD) part of the curvature of the gauge connection 

^ap = 2^aPysFap (2-3) 

(x+ =X4+xi,x_ =X4— xi) are light-cone coordinates, = p4+p\,p^ = p4 — pi) are light-cone 
momenta. 

Aw is the renormalization group invariant scale in the Wilsonian scheme defined below. 
In four dimensional Euclidean space-time, with coordinates {z,z,u,u), before the analytic con- 
tinuation to Minkowski, pi^ is the density, in units of A|,, 

p, = £5P)(z-Zp) (2.4) 
p 

of surface operators carrying at each lattice point, p, magnetic charge k and holonomy valued in 
the center, Z^, of the gauge group, i.e. such that: 

€^'^"=6'— (2.5) 

with: 

M ~ \{F^,-iF,;,)=Y^^p5^^\z-Zp{u,u)) (2.6) 
^ p 

and Zp{u,u) = Zp- The dimensionless inverse "string tension" in units of Ky^ is: 

, 10 

«' = {1.1) 
in 

with the integer, n>l, depending on the choice of the renormalization scheme, n can be absorbed 
in a redefinition of Aw. The peculiar support, (x+,x_,x+,x_), of the correlators arises as the 
projection with Minkowski signature on the base of a Lagrangian submanifold of the twistor space 
of (complexified) Euclidean space-time that occurs in our approach. 

The lattice "field of residues", jXp, and its continuum limit, jU(x), is dimensionless and normal- 
ized in such a way that the correlator in the Wilsonian scheme be renormalization group invariant. 

In the canonical scheme our result reads: 



< Tr^O^{p+,p^)Tr^O^{-p+,-p^) >£L 



{w) 

conn 



- f{-p+P-)Z-^{-p+p^) < Trj.0\p+,p^)Tr_4^0'^{-p+,-p^) >1 

(2.8) 
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with 



< Trj^O^{p+,p-)Tr,rO^{-p+,-p-) >[!fl 

~ a2 >2(2L-1) a4(2L-1) 

ti -0C'p+p-+kAl 



,2L-1 



, ( )4L-2 y ^1 

^{-p^p^r-Hog^^+... 



(2.9) 



The dots stand for contact terms, i.e. distributions whose inverse Fourier transform is supported at 
coinciding points. 

g(— and Z{—p^p^) are the (RG improved) momentum dependent canonical couphng 
and renormahzation factor ^ : 

and 

-, «^„3 I 1 „3 31ogZ 



with: 

1 11 



(471)2 3 

(2.12) 



where g = g\f^N is the 't Hooft canonical coupling constant. ^]°|^ is computed to all orders in the 
't Hooft Wilsonian coupling constant, g-w,hy: 

^ 1 10 „2 

|logZ ^ W3g>. ^2.13) 
5 log A l+cgi^ 

with c a scheme dependent arbitrary constant. Eq(2.12) reproduces the correct value of the first 
and second perturbative coefficients of the beta function. Indeed, since to the lowest order in the 
canonical coupling: 

dlogA [Any 3 

it follows: 



a log A '(471)2(471)2 3 '^"(471)2- 
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1 11 , 1 AO 44, c 



{Any 3 * (471)4 ^ 3 3 



{Anf 3 * (471)4 3 

The identification occurring in Eq(2.2) is based on the following argument. 

In the lavge-N limit there is a sector of the theory that is integrable at one loop [Q, ^ (and 
references therein), that is made by operators of ASD or SD type and their covariant derivatives. 
The corresponding anomalous dimensions can be computed as the eigenvalues of a Hamiltonian 
spin chain. 

Indeed, the anomalous dimensions of a number of operators can be computed explicitly solv- 
ing by the Bethe ansatz the Hamiltonian spin chain in the thermodynamic limit, that corresponds 
to operators of large length L and large naive dimension. In particular the anomalous dimensions 
of the antiferromagnetic ground states of length L (and naive dimension 2L) turn out to be of the 
form (see Eq.(27) of [0] and Eq.(5.23) of [g]): 

Z. = l-Lg^^^log(^) + 0(L°) (2.16) 

The ground state of the spin chain corresponds to the operators with the most negative anomalous 
dimension, that turn out to be all scalars constructed by certain contractions involving only the ASD 
part of the curvature |8||. 

For large L the anomalous dimensions of the composite surface operators 0^{p^,p-) [|l|] agree 
with the anomalous dimensions of the ground state [^, ^ of the Hamiltonian spin chain. 

The agreement is at one loop since anomalous dimensions are universal, i.e. scheme inde- 
pendent, only at one loop. Actually they agree also for L = 1, since in this case the anomalous 
dimension is determined by the beta function via the factor of g^. 

The aforementioned asymptotic identification suggests that all the glueballs in our spectrum 
are in fact scalar [|]], since this is so for the operators that correspond to the ground state of the 
Hamiltonian spin chain in the thermodynamic limit 

Thus, to summarize, there is an identification between long composite operators built by sur- 
face operators and long scalar local operators polynomial in the ASD curvature F~p, that arise as 
the ground state of the Hamiltonian spin chain, that characterizes the one-loop \arge-N integrable 
sector of YM. 

The identification is based on the coincidence of the anomalous dimensions and quantum 
numbers, by two completely independent computations. 

On one side the solution of the ground state of the Hamiltonian spin chain in the thermody- 
namic limit, i.e. the large-L limit, that involves solving an integral equation for the Bethe ansatz 

On the other side the computation of the anomalous dimensions of surface operators by com- 
bining of the Nicolai map with the localization of the holomorphic loop equation for cusped twistor 
Wilson loops on surface operators [|l|]. In particular the basic anomalous dimension for jXp arises 
by the Jacobian of the Nicolai map, while the scaling with L is a consequence of the localization. 

The two computations agree in the large-L limit. 
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Now we test our result for the glueball propagators in the large L limit against Migdal cri- 
terium. 

Firstly, we employ the perturbative one-loop RG improved expression for Zql in Eq(2.8): 

p2 y(2L) 

ZoLip") ~[iog(i^)]W 

In evaluating Eq(2.8) we employ the Euler formula, Eq(l.lO), for the spectral sum in Eq(2.9). The 
leading term, obtained substituting to the spectral sum the integral, reproduces the perturbation 
theory, Eq.(1.7), as already remarked in ^ (see last line in Eq.(2.9)). 

For the subleading contribution we obtain terms that contain the same overall Z'^^ip^) factor 
as in Eq(2.8) multiplied by inverse powers of p^, arising by the derivatives acting in Eq(2.9). 
These terms are interpreted as arising by operators of higher naive dimension, but with the same 
anomalous dimension, Z^^(/7^), as the ground state of the Hamiltonian spin chain. 

The existence of such operators is highly non trivial. In fact in [^] it has been shown that 
operators constructed by applying n covariant derivatives to the ground state have the anomalous 
dimensions (see Eq.(6.33) of |j8|]): 

7(«,L) = -7(L)+o(^) (2.18) 

Therefore in the large-L limit they are degenerate with the ground state for any finite n. QED 

Secondly, we would like to exploit the definition of the glueball propagators carrying out the 
meromorphization idea to its last consequences. From the point of view of localization of the 
cusped holomorphic loop equation [|l|] on sectors labelled by the magnetic charge, k, it is perfectly 
consistent to choose the multiplicative renormalization that occurs in the canonical scheme at the 
scale determined by the density of surface operators, pj^, i.e. on shell from the spectral point of 
view (see Eq(1.6)). 

This would match Migdal idea of meromorphization and would agree with the remark that the 
only physical information of the glueball propagator is contained in the spectral density Eq.(1.5). 
Therefore in this case, in place of Eq.(2.8), we get: 

< Tr,j.O'^{p+,p-)Tr,j^O'^{-p+,-p^) >P 



'conn 



o= 47-^ a2 7.2(2L-1) a4(2L-1) 

b\ -a'p+p-+kAl, 

" A^^z;"^ {{kAl + a'p+p_){kKl, - + {-a'p+p^ff^-' 

h\ -a'p+p-+kAl^ 

8L-8 



4^_2 y Ks\k)Z-^{k) 
h -cc'p+p- + kA^ 



with 



w 

(2.19) 



[log(-)]* 
c 
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- 1 10 
^ ~(47r)2T 

gl ~[log(-)]"' (2.20) 
c 

and the constant c to be determined by finite parts in the scheme of Remarkably the same 
leading perturbative result as in Eq.(2.8)-Eq.(2.9) is obtained: 

8L-8 

{-p^p_f^-^ f ^^^'^^ ' ^ + ... 

/t» k o (4L-4)y 

[log(-)] -'-^ Jlog(-a>+p_ + ^A^,) + ... 

- i-p^P-r-' [log(^^)] -^-"^'^^ + ... (2.21) 

But now in the Euler formula for Eq.(2. 19), as opposed to Eq.(2.8), there are not factors logarithmic 
in momentum in front of the negative powers, —j, of in the spectral sum. Thus these last terms 
cannot be anymore related, as it has been done before, to the higher dimension condensates, since 
they lack of anomalous dimensions. 

Yet, the way out now is that there are overall factors of (— )^^^^ in front of the spectral 
sum. Thus for every finite j there is always a sufficiently large L for which the product of the 
aforementioned factor and of the negative powers is indeed a positive power of p^. But these 
positive powers are irrelevant contact terms. QED 

We would like to remark that within this interpretation of meromorphization it is sufficient 
that the domain of definition of the RG flow of the coupling constant and of the anomalous di- 
mensions contains the support of the spectral measure, Eq(1.5), without the need of extending the 
RG flow to zero momentum as in [||]. Indeed a basic result of [jl], ^ is that at large-A^ there is a 
Wilsonian scheme for which the Wilsonian beta function in the ASD variables is one-loop exact 
and a canonical scheme in which the beta function has a NSVZ form that reproduces the first two 
universal perturbative coefficients Eq.(2. 15). 

This is typical of localization ^ and because of the same reason occurs also in the ^ = 1 



SUSY YM case [B^, 10 1, since in both cases the Wilsonian beta function is one-loop exact, and thus 



there is a Landau pole in the infrared for the Wilsonian coupling and, as a consequence, a finite 
range of the flow for the canonical coupling. 

But this is not a problem in Eq.(2.19) as opposed to Eq(2.8), since the RG flow needs to be 
defined only on the support of the spectral sum if we assume meromorphization. 

Finally, we discuss the interpretation, according to Migdal criterium, of correlators of surface 
operators for finite L. These are purely non-perturbative objects, because in general their anomalous 
dimensions do not match anomalous dimensions of any local operator 

Now, surface operators saturate at non-perturbative level ^ the v.e.v. of twistor Wilson loops, 
that have no cusp anomaly for large A'^. It is known that for ordinary Wilson loops cusp anomaly 
can be related to the anomalous dimension of twist-two operators (see [^] and references therein). 
Therefore the absence of cusp anomaly of twistor Wilson loops can be related to vanishing of 
anomalous dimensions of certain operators occurring in the definition of twistor Wilson loops. In 
fact the v.e.v. of these operators may actually vanish in perturbation theory because of \aige-N 
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triviality of twistor Wilson loops, due to extensive cancellations to every order of perturbation 
theory [|l|]. But these vanishing operators with no anomalous dimension may get non-perturbative 
contributions from surface operators. 

Basically the non-perturbative operators with no anomalous dimension associated to surface 
operators are powers of their lattice density, p^., defined in Eq.(2.4). 

It is only at large L that they decouple from the OPE implied by Eq.(2. 19) to match the OPE 
of the local operators that occur in the ground state of the Hamiltonian spin chain in the thermody- 
namic limit. Nevertheless, as a matter of fact, the glueball spectrum that is determined by the poles 
of the correlators of composite surface operators is L-independent and stable in the large-L Umit. 
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